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ABSTRACT 

Forms of extra compressibility terms that result from use of Favre averaging of the 
turbulence transport equations for kinetic energy and dissipation are derived. These forms 
introduce three new modeling constants, a polytropic coefficient that defines the mtere- 
lationships of the pressure, density, and enthalpy fluctuations and two constants in the 
dissipation equation that account for the non-zero pressure-dilitation and mean pressure 
gradients. 


NOMENCLATURE 

a = mean velocity of sound 

c e = modeling coefficient , Reynolds heat flux , Eq. (13) [0.35] 
( f =• modeling coefficient , eddy viscosity [0.09j 
h - enthalpy 

H — total enthalpy , Eq. (5) or Eq. (9) 
k --- turbulence kinetic energy 

k" = («)/ 2 

k = pu''u”/(2p) 

M e — free - - stream Mach number 
n — polytropic coef ficient , Eq. (10) 
p — pressure 

q 3 - molecular heat flux in the jth direction, -{ft/ Pr)hj 

Q t = Reynolds heat flux in ith direction 

Sjj = strain rate tensor , 1/2 (ujj -(- Uj t , - (2/3 )&ijU-k,k) 

T w — surface temperature 

Toe - free - stream stagnation temperature 
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T tJ = total stress tensor , viscous plus Reynolds stress 
Uj = velocity in jth direction 

(3 = fluctuating density variance , yff^ 2 /p 
Sij = Kronecker delta 
e = dissipation rate of kinetic energy 
p = molecular viscosity 
p — density 

T{j ~ viscous stress , 2fiS lJ 

A — ensemble mean variable 
A ' = fluctuating variable 

A = mass weighted ensemble mean variable 
A " = fluctuating variable in mass weighted expansion 
A yi = partial derivative with respect to time 
A >2 = partial derivative with respect to the ith coordinate 

< A > = 


INTRODUCTION 

Favre, or mass- weighted, averaging, Ref. 1, has become the most popular approach 
for establishing the field equations used in evaluating compressible turbulent flows. One 
reason for this is that the equations for mass, momentum and energy conservation retain 
the relative simplicity, term for term, of their incompressible counterparts, Ref 2. The ef- 
fects of density fluctuations are submerged within the definitions of the Reynolds stresses 
and heat or mass fluxes and the profusion of new terms that occurs when primitive vari- 
ables are employed, e.g. Ref. 3, is avoided. The effects of density fluctuations, however, 
become reintroduced into the problem when use is made of field equations representing 
the turbulence moments and scale. These equations contain additional terms unique to 
compressible flows that involve density fluctuations and the non-zero divergence of the 
turbulent fluctuating velocities. The purpose of this paper is to reexamine the modeling 
these compressibility terms, with emphasis being given to the terms that appear in two- 
equation models of turbulence. Of course, similar forms of these compressibility terms can 

be applied to second-order closure models as well. 

The possible importance of these compressibility terms was first suggested in the paper 
by Wilcox and Alber, Ref. 4. These authors modeled the aggregate of these terms in the 
turbulence kinetic energy equation with an analysis that employed the assumptions that 
pressure and total temperature fluctuations can be neglected. They developed a model 
for these terms, expressed in mean variables, that introduced a single additional modeling 
coefficient. No corrections were applied to the scale equation. The turbulence model, with 
a particular modeling coefficient, was applied to some selected supersonic flows, including 
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near wake boundary layers with Mach numbers up to 8. The computed results agreed 
resonabley well with the experimental data, but the absence of computed results without 
the additional compressibility term in this reference did not permit firm conclusions to be 
drawn by a reader as to t he value of including this term. 

Oh, Ref. 5, reexamined these compressibility terms in the turbulence kinetic energy 
equation with the view of improving computations of the spread rate of a free-shear layer 
formed by a single supersonic stream. As in Ref. 4, no modifications were made to the 
scale of turbulence, which, in this case, was prescribed algebraically. Oh s analysis of 
the compressibility effects was based on an assumed detailed structure of two-dimensional 
shocklets developed by the turbulent eddies in relative supersonic motion. In addition, a 
form of Taylor 's hypothesis for the turbulent motions was invoked. Although, the process 
introduced two new modeling coefficients, one could be neglected by Oh as it was a factor 
of the term involving the mean pressure gradients, which were set to zero in the free-shear 
problem. The second modeling coefficient was shown to produce a diminished spread 
rate with increased Mach number, as occurs within the experimental data. It was found, 
however, that good agreement with the data over the range of Mach number from 1 to 
4 would have required altering the value of the coefficient by about 40% over the Mach 
number range. 

To avoid conjectures, unsupported by experimental data, regarding the detailed struc- 
ture of shock waves created by colliding three-dimensional eddies, Rubesin, Ref. 6, derived 
formulations for the extra compressibility terms in the turbulence kinetic energy equa- 
tion by introducing the concept that the state conditions within an eddy could be related 
polytropically. This was another way of accounting for non-isentropic processes, shocklets, 
etc., that occur within the turbulence, but without introducing structural assumptions. It 
was necessary, however, to make assumptions regarding the fluctuation of the total tem- 
perature and the variance of the density fluctuations. In Ref. 6, it was assumed that the 
total temperature did not fluctuate and that the variance of the density fluctuations varied 
very slowly along streamlines. The resulting formulations for the compressibility terms in 
the kinetic energy equation contained a single modeling coefficient, namely, the polytropic 
coefficient . 

The work of Ref. 6, was extended by Vandromme, Ref. 7, to retain the effects of 
large mean density variations across streamlines and to introduce extra compressibility 
terms into the turbulence scale equation, i.e. the rate of dissipation of turbulence kinetic 
energy, as well as in the kinetic energy equation. In addition to the polytropic coefficient, 
Vandromme introduced a modeling coefficient in the added term in the dissipation rate 
equation. Vandromme found that the spread rate of the free-shear layer in supersonic flow 
could be represented very well over the Mach number range from 0 to 4 with single values 
of the polytropic coefficient and the modeling coefficient in the dissipation equation. 

There have been several papers addressed to the evaluation of compressibility effects 
within internal combustion engines, e.g. Refs. 8-10. In Refs. 8 and 9, the burden of 
maintaining reasonable length scale behavior with a two-equation model during volumetric 
compressions was placed on the dissipation equation, with the result that an an additional 
divergence term was added to the dissipation equation in k-t turbulence modeling. The 
turbulence kinetic energy equation was left without additional terms, which, as argued in 

3 

ORIGINAL PAGE IS 
OF POOR QUALITY 



Ref. 10, is consistent with the low Mach numbers characteristic of internal combustion 
engines. The difference, then, within this group of papers is in the values chosen for the 
modeling coefficient multiplying the divergence term. 

In a recent paper. Ref. 11. Zeman identified the u dilitation dissipation as the 
extra compressibility term in the kinetic energy equation. In a manner similar to Oh, 
Zeman utilized a modeled shock wave structure within the turbulent eddies to obtain the 
functional form of the compressibility correction. Assumptions regarding the shock wave 
thickness relative to the eddies and the probability distribution function of the turbulence 
Mach number led to a form of the compressibility term that contained two new modeling 
coefficient s. one of them being t he kurtosis of t he t urbulence. It was shown that appropriate 
constant values of these modeling coefficients in computations of the behavior of a free- 
shear layer could be made to represent the data of Ref. 12 quite well. 

From this brief review, it is apparent that various authors chose considerably different 
routes to introduce compressibility terms into the field equations for either the turbulent 
kinetic energy equation or the dissipation equation, or both. Under high Mach number 
conditions, the order of magnitude arguments in Ref. 10 would suggest that these terms 
should be included in both equations. 

Because of the flexibility introduced into these models by the introduction of one 
or more additional modeling coefficients, agreement with particular experimental data, for 
example, high-speed free-shear experiments, could be achieved by several disparate models, 
Refs. 5, 7, and 11. In the absence of the application of these same models to more complex 
flow situations, such as flows where shock waves impinge on boundary layers to produce 
separation or where they interact with shear layers away from surfaces, it is not known 
if the various models will differentiate and show one that is superior. To perform such a 
comparative study, it will be best to utilize a single computer code that contains a series 
of the alternative turbulence models cited here which can be turned on and off at will. 

The present paper is devoted to the development of one such model. It is an extension 
of the work of Refs. (6) and (7), but with particular attention being given to the non- 
adiabatic character of high- Mach number flow fields and their large variations in the fluid 
properties. In addition, the model allows inhomogenieties in the turbulence. 

ANALYSIS 

When Favre or mass- weighted ensemble averaging is applied to the conservation equa- 
tions describing the turbulent flow of a compressible fluid, the following set of equations 
results. The only physical assumptions applied to these equations are that the fluctua- 
tions in the fluid viscosity and thermal conductivity can be ignored. Turbulence effects 
on molecular properties are not expected to be important where the molecular properties 
themselves dominate. 

Continuity 


p,t + (pwjKj = 0 


(i) 


Momentum in the ith direction 


4 



( 2 ) 


(pUr) y t -+ {pUjU i + SijP - Tij)j - 0 

Here, the total stress is composed of molecular and turbulence contributions. 


T 


1 2J 


4 - 


T rj ~ P U i U j 


(3) 


Total energy 


(pH - p), t + (pUjH)j ~ 


\pu 


I’jh" + pu'jk" + qj + q'j - UiTij ~ u\'T t 


u”t" 

l IJ 


(4) 


The quantity H represents the mean total enthalpy, which is the sum of the mean 
static enthalpy, h, and the kinetic energies of the mean motion plus the turbulence. 

H = h + UiUif 2 + pu"u" / (2 p) (5) 

In addition, the turbulence kinetic energy equation can be written as 


(pk) t + (pujk)j — 


- (pu'jk" + t I'jp - <t'' )j 


-t Wj + + p'u'r - u'fpj 


( 6 ) 


The additional compressibility terms in these equations, absent for incompressible 

fluids, are those involving u", r", q" and u",. 

Before evaluating the quantities that contain these additional compressibility terms, 
it is worthwhile to eliminate those of little consequence. In regions where the molecular 
transport properties are important, it is to be expected that the_effects of density Aqu- 
ations should be insignificant. With this assumption, the terms r,", q " , v"f i} , and v,fij,j 
can be omitted from Eqs. (3), (4), and (6). 

These equations then reduce to 


T,j = r,j - pujuj 


(pH - p), t -r ( pUjH = 


pu'jh" + pu'jk" + qj u i Tj j u ■ r. 


n „tt 

ij 


and 


Q?d2:y 
OF FC 


L L h. 




<*• 'i/i > 



( 9 ) 


(pk) t + {pujk)j = - (pu"u'j)ui,j ~ (pu "k" -r u'jp 

+ 7<, - K'P.i 


u"t" 


)J 


The first step required in the modeling of the remaining additional terms is to deter- 
mine a relationship between the fluctuations of velocities and the fluctuations of the state 
properties of the fluid. This is initiated by following Ref. 6 and assuming that the fluid 
behaves in a polytropic manner. Thus, the polytropic coefficient, n, becomes a turbulence 
modeling parameter. 

Expressed in mass weighted variables, this assumption becomes 

p' n p' _ ” pT _ n ph" 

p p n — 1 pf n — 1 ph 

The last relationship is only exact for constant values of the specific heat, but will be 
used here as a first approximation even when the specific heat does vary. 

Next, a relationship between T" or h" and u" has to be established. In Ref. 6, it was 
assumed that the total enthalpy at a point in space remained constant with time despite 
the turbulence. This is expressed as 

h” = h H + u i u , ; + (^-~k) = o (ii) 

If the higher order terms are neglected, then 

h" = — Uiu" (12) 

The experimental results cited in Ref. 13, however, indicate that within an hypersonic 
boundary layer, e.g. at M e - 6.4 over a cooled surface where T w /T 0e = 0.46, the local 
r.m.s. total temperature fluctuations are not zero, but can have magnitudes ranging from 

1 .4% < < 4.8% 

H 

with the larger variances occurring closer to the surface. Equation (12), therefore, has to 
be modified to account for fluctuations in the total temperature. 

One could expect that the local temperature or enthalpy fluctuation would be de- 
pendent on the intensity of the turbulence and, also, on the mean static temperature or 
enthalpy gradients. One such form, analagous to a mixing length formulation, is 

h" = —au"h t i (Id) 

where a is a local mean parameter, yet to be determined. 

If Eq. (13) is multiplied by pit" and averaged, there results 

Qj - pu"h" = -apu"u"h t j (14) 
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which relates the local heat flux to the Reynolds stresses and the mean enthalpy gradients. 
This form is consistent with the modeled forms developed in Ref. 14, namely 


with 


c c = 0.35 


Consistency between Eqs. (14) and (15) requires 


a -■ 


c e k/i 


(16) 


so that Eq.(13) becomes 



A more tensorially complex form of Eq. (13) could have been adopted to be consistent 
with the more recent, and empirically justified forms of Ref. 15. However, since the current 
work of introducing the effects of compressibility into the k and f equations will later involve 
thin layer approximations which force the principal fluxes of Refs. 14 and 15 to be identical, 

the simpler form of Eq. (13) will be used here. 

The experimental data of Ref. 13 can be used to test, the validity of Eq. (17) provided 
the differences in the physical assumptions and dependent, variables used ill the data reduc- 
tion and this analysis can be reconciled. The data are expressed in unweighted ensemble 
averaged quantities whereas mass weighted variables are used here. The quantities mea- 
sured in the experiment through the use of a hot wire with a large range of overheat values 
were the fluctuations of total temperature and axial mass flow. These quantities were then 
converted to static temperature and velocity fluctuations through the assumptions m Ref. 
13 that static pressure fluctuations and higher order correlations could be neglected. With 
the latter assumption, the distinction between mass weighted and unweighted ensemble 
averaged quantities vanishes. Consequently, a direct comparison between the results of 
Ref. 13 and Eq. (17) is possible. 

In a boundary layer, the principal terms on the right side of Eq. (1 7) are 


h" = 



(18) 


Similarly, the principal shear stress in the boundary layer, given by two equation modeling, 

r i2 = C j l)ui 

When Eq. (19) is used to eliminate c in Eq. (18), the variances of the static enthalpy 
and the velocity can be related as 


< h" > jh 

< v" > ,/fi] 



( Ti2 / p) <U ft 
k <«" 


2 > 
1> 


(fil/wif) 


d(h/h r ) 
d{ it i /u, f ) 


(h/h c ) 


( 20 ) 



where the differentiation is performed at a fixed station. 

In evaluating the right member of Eq. (20). the following values are employed: 

c t - 0.35 
C\ = 0.09 

( J- )/k = 0.3 

P 

which is consistent with the above value of C\, and the universal log law', and 

< u'.l > / < u'l >- 0.64 

which corresponds to the simulated data of Ref. 16 at a y of about 140. The mean 
enthalpy and velocity relationships are taken directly from the data of Ref. 13. The 
results from Eq. (20). winch is based on Eq. ( 1 t ), are compared w'lth the corresponding 
experimental data of Ref. 13 in Fig. 1. In addition. Fig. 1 show's the relationship of the 
enthalpy, velocity variances that result from the use of Eq. (12) as evaluated in Ref. 6. 

The use of Eq. (17) results in a ratio of the enthalpy and velocity variances that 
represents the hypersonic boundary layer experimental data of Ref. (13) reasonably well 
over most of the boundary layer. It is surprising that the computed results differ from the 
data as much as they do in the log region of the boundary layer, y/6 < 0.2, where the 
assumptions of the anisotropies and shear stress relationships listed above whould have 
been expected to apply the best. In the region near the boundary layer edge, how'ever, it 
is not surprising that the computed results are high because no allowance has been made 
for the shear stress to approach a zero value there. Nonetheless, these results based on 
Eq. (17) are a significant improvement over those that result from the use of Eq. (12), 
the basis of the w'ork in Ref. 6. Accordingly, Eq. (17) is used in subsequent steps of this 
analysis. Note that Eq. (17) indicates that the sign of the effect on the enthalpy fluctuation 
caused by a velocity fluctuation depends on the sign of the local mean enthalpy gradient. 
Alternative directions of heat flux appear to have opposite effects on the turbulence kinetic 
energy. 


Evaluation of u ’■ 


From Eqs. (10) and (17 


P_ 

P 


1 ph" 

n 1 ) ph 


1 P k h , , „ 


( 77 — 1 ) ' p t h 3 
which can be used with the definition of v'' to yield 


P' u " 


(n - 1) p f h 


- e k h j - 


( 21 ) 


( 22 ) 
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Recall that 


(23) 


h = 


a 2 


(7-1) 


so that Eq. (22) can also be expressed as 




pa 2 


(24) 


The bracketed terms represent moments of turbulence Mach numbers, which vanish 
in incompressible flow, i. e. a -> oo. Thus u" is a measure of the degree of compressibility 
of the turbulence. It is also interesting that i/'' has a value of zero in the absence of any 
local heat transfer, where hj = 0. The latter point can be seen directly when Eqs. (24) 
and (15) are combined to yield 


(7-i) Q> 
( n — 1 ) pci 2 


(25) 


In most aerodynamic applications, the principal mean pressure gradient is p 1? which 

identifies u" as the main quantity required in the last term on the right of the kinetic 
energy equation, Eq. (9). If the flow also behaves as a thin shear layer where h , 2 is much 

greater than h tl or h, 3 , Eq. (24) reduces to 


(7 


1 ) 

~ c ( , — h 7 
72 - 1) ‘ f 


P u " u 2 

2 

pa* 


(26) 


Since the Reynolds sliear_stress and the static enthalpy gradient in Eq. (26) can 
each possess different sign s, u" can be either positive or negative. For example, within 
a planar, hot jet, h )2 and P u " v 2 c hange sign simultaneously on the centerline of the jet. 
Thus u[' retains a negative sign on either side of the centerline. Alternatively, in a high 
speed boundary layer over a cooled surface. ~h has a maximum value within the boundary 
layer, while the shear stress maintains the same sign everywhere. Here, the sign of v" is 
different below or above the point wdiere h is maximum and the k equation will be affected 
differently in these different zones of the boundary layer. This may explain, in part, why 
attached boundary layers seem to exhibit less compressibility effects than do free shear 
layers at similar Mach numbers. 


Evaluation of p'u" ■ in the Kinetic Energy Equation 


In terms of mass weighted variables, the continuity equation for the fluctuating quan- 
tities is 

Pp - (P'v-j + P u 'j),j = 0 ( 2 ' ) 
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After some algebraic manipulations and neglect of higher order terms, Eq. (2/) becomes 
(p' 2 ),t - Uj{p t2 ),j 4- 2 p l2 Ujj -4 2 PjP'u'j -r 2pp’u" J = 0 (28) 


or 




'jo 




i r 


(p' 2 ), t 4 - 


P p J 2p 

Finally, relating the pressure and density fluctuations as in Eq. (K), allows writing 

r><> = j^h 

or, with Eq. (28), the pressure rate of strain term becomes 

[(p' 2 ),< - »,j(p' 2 ),j 


(29) 


(30) 


-nr -1 - p’ 2 
P V i,j = n P\~ U ^J 2 


P P 2 p 2 


(31) 


Eq. (31) shows that the variance of the density fluct uations, ,3 = yj p 12 /p, plays a role in 
establishing the magnitude of the pressure, rate of strain term appropriate to the kinetic 
energy equation. In fact, when Eq. (31) is expressed directly in terms of_tdie variance 
0, and use is made of the mean continuity equation and the definition of u " , Eq. (22) 
simplifies to give 


pV' i = np 


- J P’ J 


1 


}<) 


U * 2 


(3 2 ),t + «j(/7 2 )o 


(32) 


Locally, the density variance can be established from Eqs. (10) and (17) as 

3 2 ■ 


V ( " - 1 )ph 

or with higher order terms dropped 


— | puVvfh/hj] 
- \ )ph « ) 2 


(33) 


.P 


c t . 


k\ 1 


(n — 1 )h 


ikmmj 


(34) 


In a thin shear layer, where h <2 predominates, Eq. (34) reduces to 


3 2 


n o 

P»2~ 


1 ) 2 h 2 f 2 \ pk 



(35) 


In the computation of a typical k-t model, each of the quantities appearing in Eq. (35) 
are known local quantities and can be used to establish the local field value of the density 
variance required in Eq. (32). 
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The remaining terms in Eq. (9) that require modeling are the transport or diffusion 

of kinetic energy 

(36) 


D k = -(pu'-k" + u'jp' - 


and the turbulence dissipation 




(37) 


The reasons for identifying the whole of Eq. (36) as a transport term is its divergence form 
and its appearance in the total enthalpy equation, Eq. (8), where the pressure fluctuation is 
an implicit part of the enthalpy fluctuation. The dissipation term, Eq. (3i ), however, does 
not occur in the total enthalpy equation since it merely represents an exchange between 
mechanical energy and heat, the sum of which contributes to the dependent variable of 
that equation. In incompressible flows, and in Ref. (11), these terms are combined and 
their difference is modeled. Here, however, these terms will be treated separately because 
their molecular terms are important in different parts of the boundary layer and this can 
be used to advantage in the turbulence modeling of the additional compressibility terms. 

The turbulent parts of D k are modeled to be consistent with the level of modeling in 

Eq. (15), namely 

(38) 


pu"k" + u"p = -c.-pu'ju'jk'i 


where c s = 0.25, Ref. 17. . 

For a thin shear layer, the principal direction of gradients m the kinetic energy require 

i = j = 2, so that Eq. (38) reduces to 


pu”k H + u'V 


Cs-Pu"u"k ' 2 

€ 


(39) 


With C S = 0.25 and pu'^/pk ~ 0.36, Eq. (29) can be rewritten as 


p 

pu'^k" + u'V - -Mpjkp (40) 

which is consistent with the simple gradient diffusion form of the standard k - e model 
formulation. 

The viscous part of Eq. (36) is 


where the effects of the fluctuating part of the viscosity have been neglected. Because 
these viscous terms are only important in comparison with their turbulent counterparts 
verv close to the surface where the turbulence has been damped, the extra compressibility 
effects resulting from mass weighting are not expected to be important where the viscous 
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terms dominate. The turbulence, therefore, in these terms can be treated as solenoidal 
and u'V//2 - fc. Eq. (41) then becomes 




(42) 


In incompressible flows, when the diffusion and dissipation terms are combined, only the 
first term on the right of Eq. (42) remains to define the viscous diffusion of the turbulence 
kinetic energy. This turns out to be a reasonable approximation even in inhomogeneous 
compressible turbulent flows. In the near wall region, the principal gradients of the tur- 
bulence moments are normal to the surface, where i = J - 2. At the surface, k : ~ y , 
whereas uW ~ y\ so that the latter can be neglected. The second term becomes a larger 
contributor with distance from the surface, however, there the contributions of both terms 
in Eq. (42) become less significant in comparison with the turbulent transport. It has 
been estimated with simple mixing length arguments that the second term in Eq. (42) 
increases the local transport of kinetic energy by a maximum of 13 percent at a value of 
y - - 12. Over the entire sublayer, 0 < y * < 60. this term increases the diffusion of k only 
by about 4 percent , so that Eq. (42) can be simplified to 





(43) 


which can be used in both Eqs. (8) and (9). 
With Eqs. (38) and (43) 


lh 


* j 


(44) 


or in a thin layer 



( 


(45) 


The dissipation term, Eq. (37), expands to 








y» 


X, 




j>K. 


(46) 


It is noted that the “dilitation dissipation 5 ', defined in Hef. 11 and modeled there, 
appears as the third term on the right of Eq. (46). Here, to account for inhomogeneous 
turbulence, the dissipation will be treated as consisting of all three terms in Eq. (46) and 
its difference from the first term, the < of the standard k — t model, will require modification 
of the form of the e equation. 
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Those readers familiar with the details of Ref. 11 will note that the “dihtation dissi- 
pation” term above has a different sign and coefficient from that of Ref. 11. The reason 
for this is that the form of the total dissipation of Ref. 11 resulted from manipulation of 
the second term of Eq. (46) under the assumption that the fluctuating turbulent velocity 
moments, not necessarily the Reynolds stresses which include local density, are homoge- 
neous. It is interesting that an identical result occurs if the diffusion and dissipation terms 
are first combined, but then it is required to assume that both the velocity moments and 
the viscosity are homogeneous. The approach used in this paper avoids these restrictions, 
and allows for both inhomogeneous turbulence and physical properties that are variable. 

In summary, the k - e model has been modified for the additional effects of com- 
pressibility in the following manner. The kinetic energy equation, Eq. (9), is rewritten 

as — _ 

(pk) t + {pujk) tj = -{puWjfaj - p( + D h 4 p u'l, - (47) 

where D k is given by Eq. ( 44 ), pVE is given by Eqs. (32) and (34) or (35), and w" is 
given by Eq. (24) or Eq. (26). The dissipation, e, is modeled with the field equation 


(/«),« + 




-0 



K 


(48) 


where 


k 




(49) 


To most clearly illustrate the additional terms introduced to account for the effects 
of compressibility, equations (47) and (48) are given in their high Reynolds number forms, 
absent any near wall corrections. Eq. (48) contains three terms in addition to those 
contained in the c equation of the standard k-e model. The third term on the right has 
been added to account for the expected dependence of turbulence length scale on passing 
through a shock wave following arguments similar to those of Refs. (9); it does not contain 
the factor of 1/3 recommended in Ref. (10) because the latter is only appropriate to 
isotropic turbulence. This term is necessary even when the effects of fluctuations in density 
are ignored. The other two terms, those containing C (3 and C f4 , account for the effects of 
the fluctuating density on altering the total dissipation of kinetic energy. In addition to 
the polytropic coefficient , n, the two new modeling constants introduced are the ( , 3 and 
C (4 . All the other modeling constants and relationships between k , c, the eddy viscosity, 
etc. are the same as in the usual incompressible k-e approaches. 


CONCLUDING REMARKS 


The forgoing anlysis to define the extra compressibility terms that are introduced 
by fluctuating density in the Favre averaging process has resulted in the introduction of 
three new modeling constants: the poytropic coefficient, n, and two coefficients in the 
modified kinetic energy dissipation equation, (\ 3 and C\ 4 . In the absence of validation 
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computations that include these terms for experimental conditions where they can be 
important, it isn’t possible at this time to give recommendations for their appropriate 
values. In Ref. 6 , it was recommended that n = 1.2. The basis for this recommendation 
was no more profound than that n = 1 .2 is consistent with polytropic irreversable expansion 
relationships in thermodynamic cycles often cited in standard thermodynamic text books,, 
e.g. Ref. 18. On the other hand, when the behavior accross a normal shock wave is 
considered to be polvtropic. it is found that n > q. the isentropic coefficient. In view of 
these differences and in the absence of any empirical guidance, it is believed that starting 
with n— 1.4 is reasonable. It is the intent of the author to establish the value of C f 3 first, 
through comparisons of computations with experimental data for flow fields involving zero 
mean pressure gradients, thereby eliminating the term containing C * 4 . C andidate flows for 
this are single or double stream free- shear layers and attached flat plate boundary layers, 
all in high speed flows. The effects of the added terms are expected to be less for the 
attached flow. Once consistent values of n and C f e 3 are established, test cases involving 
strong pressure gradients will be pursued to establish C f 4 . 
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